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We propose to apply "worldline numerics" to a numerical calculation of quark determinants. The 
Gross-Neveu model with a U(l) chiral symmetry is considered as a first test. The worldline ap- 
proach allows for an analytic renormalisation, and only finite parts of the determinant require a 
numerical calculation. It is shown that the discretisation of the worldlines, which is central to the 
numerical treatment, preserves chiral symmetry exactly. Numerical results for a kink configura- 
tion as a scalar background field are shown and compared with analytical results. The case of 
finite fermion chemical potential is also briefly discussed. 
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1. Introduction: 

Over the last decade, the QCD phase diagram as a function of temperature and baryon density 
has attracted intense investigations using computer simulations and collider experiments such as 
RHIC undertaken at the Brookhaven National Laboratory. From the simulation point of view, 
many efforts in lattice gauge theory in the recent past were devoted to control the severe sign 
problem. The proposals use a Taylor expansion with respect to the baryon chemical potential 
^ ||], imaginary values for the chemical potential /i [Q, or overlap enhancing techniques |^]. 
Despite of these successes, our knowledge of the QCD phase diagram is still limited to rather small 
values of the chemical potential. 

Further insights have come from QCD-like theories such as 2-colour QCD [Q, ^, |T^], from 
perturbative studies and QCD motivated models. In the latter case, one assumes that the phase 
at highest densities is homogeneous and argues that quark matter is organised in a colour super- 



conducting state Qlll , |12| , |13| , |14| ]. Recent studies of the Gross Neveu model in the limit of many 



flavours have attracted a lot of interest since it was found that the high-density state of fermion 



matter forms an inhomogeneous "baryon crystal" UlSQ. 

Good chiral properties of the fermion action is of central importance for an investigation of 
quark matter at intermediate densities, since the high-density transition is driven by chiral dynam- 
ics. Unfortunately, lattice fermion actions necessarily suffer from the fermion doubling problem as 
firstly pointed out by Nielsen and Ninomiya [|^]. Nowadays, staggered fermions [[T7|], domain wall 



fermions [[1^] or Neuberger fermions [[191], which are an explicit realisation of the Ginsparg-Wilson 
relation [20], are widely used in numerical simulations. Despite of these advanced formulations 
and great numerical efforts, it turns out cumbersome to achieve good chiral properties such as a 
sufficiently small pion mass. 

Since the worldline approach to the quark determinant does not use a lattice discretisation 
of space-time, it circumvents many of these significant difficulties. Here, we will argue that the 
prospects of the worldline approach are (i) exact chiral symmetry but yet a fully numerical ap- 
proach, (ii) analytic renormalisation and (iii) a clear description of Fermi surface effects. 



The worldline method is a string-inspired approach to quantum field theory; see [ ]21[ ] for a 
review. It was further developed into a viable tool for an efficient calculation of functional de- 
terminants for arbitrary background fields [p2|]. Subsequently, worldline numerics has enjoyed a 



wide span of applications ranging from the Casimir effect [ |23| , [24[ ] and fermion induced quantum 
interactions [ p5[ ] to the description of pair production in inhomogeneous fields [26|. A worldline 
lattice formulation has been presented in [p7|]. 

2. The chiral Gross-Neveu model 

2.1 Setup of the model 

The Gross-Neveu model in its original formulation is a two dimensional fermionic theory 
which shares with QCD the property of spontaneous chiral symmetry breaking and asymptotic 



freedom [28|. Due to Thies and Urlichs, the phase diagram is analytically known in the limit of 



many flavours [15|. This model therefore provides for a benchmark test for any new numerical 



method which tries to extend its reach to very dense fermionic systems. 
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In the chiral version of this model, a pseudo-scalar field k{x) acts as chiral partner of the scalar 
field o{x). The partition function is given by 



^ = j SO&K exp{-A^5fer - A^5bos} , 5bos = ^ j [^^W + 



-trln 
2 



, (2.1) 
(2.2) 



where g is the bare coupling constant and where we have used anti-hermitian Dirac matrices. The 
partition function is invariant under a ?7(1) chiral rotation of the fields: 



a'(x) 



COS0 — sin0 
sin d cos d 



o{x) 
7l{x) 



(2.3) 



In the \wge-N limit, one assumes that fluctuations of the mesons are negligible, and that the relevant 
field configurations can be obtained in leading-order saddle point approximation: 



Sfer + d X 



2gn 



o^{x) + n^{x) 



mm. 



(2.4) 



2.2 The worldline approach to the GN model 

The key ingredient of the worldline calculation of the fermionic part 5fer in (|2.2[) is the repre- 



sentation of Sfer in terms of an ensemble average of closed loops (t), T = . . . T , (0) = (T), 
in Euclidean spacetime. In the loop cloud approach [22|, the worldlines are generated according to 
the free probabiUstic measure: 



5(xcm[x]-Xc^ expj- dz ^ I , 



where the loop centre of mass given by 



is constrained to Xc. The fermion determinant is then represented by 
r(a, k) = try ^exp ^/ J dx {^o + ijs ^k) ^ , 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



where A is a UV cutoff. 

2.3 Exact chiral symmetry 

In numerical calculations, a closed loop x{t) is represented by a finite number of points: 

Xi x{ti), i=l...Np, dT = T/Np . 
The spin factor r(a, n) is approximated by a path-ordered product 

rdis(a,7r) = tr^ Yl^exp{idT [^a(x,-) + iysMxi)] ) , (2.9) 
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The crucial observation is that, in spite of the discretisation, Fjis still is exactly chirally invariant. 
To show this, we define a unitary matrix U hy U = cos(0/2) + /sin(0/2) 75 , and show that 

^o'{x) + iys^n'ix) = U {^o{x) + iy5^n{x)) uK (2.10) 

Because of the path ordering and the closeness of the (discretised) loops, we easily find that (see 
figure [1], left panel for an illustration) 

^Ais{o\%') = rdis(a,7r) . 



The other parts of the fermionic action ( |2.7| ) as well as the integration measure for the mesonic 
fields are trivially invariant (there is no anomaly in this model) leaving us with an exact chiral 
symmetry for the discretised theory. 

2.4 Renormalisation 

Another big advantage of the worldline approach to fermionic determinants is that the UV 
regularisation can be performed along the lines made explicit in the ab initio continuum formula- 
tion. Only finite parts of the determinant must be calculated by numerical means. This implies that 
one does not need to invoke any "order-c?" improvement which is instrumental when conventional 
lattice fermions are considered. Let us illustrate the renormalisation procedure for the present case. 
Introducing the space-time average 

= cFx [o'^{x) + Ti^{x)], (2.11) 
the fermionic action can be split into a UV divergent and a finite part: 

5fer = 5o(M,A) + St[a,%] , (2.12) 

-2exp{-rM2} ^ , (2.13) 

r2 ^00 j'p 

5o(M,A) = - J — expl-TM^} . (2.14) 



With this construction, the part S^"^ of the action which involves time consuming numerical simula- 
tions is UV and IR finite. Accordingly, we have removed the regulator in (2.13) by taking the limit 
A ^ 00. The part Sq of the action contains the divergent pieces which can be calculated explicitly; 
dropping a field-independent constant, we obtain 

5o(M,A) = — ym'-Xn^ + {ye-\)m'-\ + ^( , (2-15) 

where Je is Euler's constant. Adding the bare bosonic part of the action in (2.1), we can impose 
renormalization conditions, for instance, of Coleman-Weinberg type; this defines the renormalized 
coupling at an RG scale ji, g^^ijJ.) := d^S/do^\j^2^(j2^^2, finally yielding, 

5o(M,A) + 5bos = ^ m2 fln^ - 1 V = (2.16) 

471 \ Mq / 
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O worldline 100 ppl 
exact: sech{x)'^2 




Figure 1: Illustration of chiral invariance of the spin factor (left); zero-mode wave function reconstructed 
from the heat kernel using free loop ensembles (right). 



where we have traded the coupling for an RG invariant mass scale Mq in the large N hmit, 
reflecting dimensional transmutation. This scale also denotes the large-A'^ minimum of the action at 
zero temperature and density, = Mq =const. 

2.5 A numerical benchmark test 

For a benchmark test, we choose a kink configuration as a background field: 

a(?,x) = a(x) = tanh(x) , n{t,x) = 0. 
This kink interpolates between the two homogeneous vacuum states a = ±1 and is the basic build- 



ing block of the 'baryonic crystal' of the Gross-Neveu model [ 15]. The Dirac structure decomposes 
into two Schrodinger problems for the heat-kernel traces 

tr,exp{-r//±} = ^^/exp{-r//±}\ , H± = -d'^ + ±^o{x). (2.17) 
yAiiT \ I X ax 

On the other hand, the heat kernel expectation value can be expressed in terms of the eigenmodes 
of the Hamiltonian H±. In particular, H = —d^ + 1 — 2/cosh^;c gives rise to a zero mode, 

-^(exp{-r//_})^ = Iv^oWl' + IIV^/Wl'e-^^'. (2.18) 

This choice therefore challenges the worldline approach, and the crucial question is whether a 
moderate number of free loops is able to grasp the zero-mode contribution. We have calculated 
the heat kernel expectation value for a range of propertime values T of ^(10) using 50000 loops 
consisting of 100 points per loop, all of which include the point x. In this propertime range, the 
contributions of the excited states is small, and the modulus of the zero-mode wavef unction Iv^ol^ 
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can be reconstructed from (2.18) by a fit. Our numerical findings for the modulus of the zero 
mode wavefunction are compared in figure |l|, right panel, with the exact result. A rather rough 
discretisation of the loops already yields quite accurate results. 



3. Finite densities 



Let us consider the case of a time independent, but non-homogeneous scalar field {k is set to 
zero for illustration purposes), which is relevant for the high density crystal phase. The fermion 
determinant at finite temperatures l/j8 and finite fermion chemical potential ji in the worldline 
approach is given by 



1 



dT 



L exp. 



271 



(?i+l/2) -ijx 



K{T), 



2 7i/A2 T 

K{T) = ^ J dx^ ^expj-^^ Jta^l tr^^cosh^/^ dx 



(3.1) 
(3.2) 



In order to make the physics of the Fermi surface transparent, we introduce the Laplace transform 
of the kernel K{T) hy 

poo 

(3.3) 



K{T) = 2 / dE E exp[-TE^] p{E) , 
Jq 

where p{E) has the interpretation of the density of states. We obtain: 



— (n+1/2) -iAt 



1 2 



TE' 



(3.4) 



The technical advantage is that we have mapped the problem of dealing numerically with the Fermi 
surface onto the problem of a free particle theory with single particle energy E. Hence, it is well 
known how to evaluate the proper time integration T and the Matsubara sum n in (^^. Decom- 
posing the fermionic action into temperature dependent and independent parts, we find: 



^temp 
•^fer 



^temp ^0 

dEEp{E) 



In 



1 + exp{-/3(£ + At)} 



+ In 



1 +exp{-j8(£-At)} 



S^^ = ldEEpiE)\[/-Lli'^ exp 



1/A2 T 



2n 



{-T [kl+E^]'}. 



(3.5) 
(3.6) 

(3.7) 



Hence, the formulation offers a complete control over the physics associated with the Fermi sur- 
face, and even the low temperature and high density regime is accessible. For instance, the small 
temperature expansion (arbitrary chemical potential) of the baryon density is given by 



1 d 



temp 



j^dEEp{E) + A [Ep{E)] \e=^ + e{T') . (3.8) 



Fermi surface effects can thus be studied in a systematic fashion. Note, however, that the numerical 



calculation of the density of states p{E) from the kernel K{T) in (3.3) can be cumbersome 
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